Faraday Wave Pattern Selection Via Multi-Frequency Forcing 
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We use symmetry considerations to investigate how damped modes affect pattern selection in 
mult i- frequency forced Faraday waves. We classify and tabulate the most important damped modes 
and determine how the corresponding resonant triad interactions depend on the forcing parameters. 
The relative phase of the forcing terms may be used to enhance or suppress the nonlinear interac- 
tions. We compare our predictions with numerical results and discuss their implications for recent 
experiments. Our results suggest how to design mult i- frequency forcing functions that favor chosen 
patterns in the lab. 

PACS numbers: 05.45.-a, 47.35. +i, 47.54.+r, 89.75.Kd 
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Faraday observed in 1831 that patterns of subharmonic 
standing waves form on the surface of a fluid when its 
supporting container is shaken vertically with increasing 
strength [l| . The first Faraday wave experiments (see |2j 
for a review) used sinusoidal forcing to shake the fluid, 
and produced simple patterns (e.g., stripes, squares, and 
hexagons). Remarkably, this experiment remains a rich 
source of intriguing new patterns, most recently as the 
basic sinusoidal forcing has been replaced by periodic 
functions with two or three frequency components. The 
two-frequency case has been investigated a great deal ex- 
perimentally 0, 0, S IE an d theoretically @, 0, 0, El , 
and has produced a variety of states, including super- 
lattice patterns and localized structures, as well as the 
first experimental realizations of quashoatterns. Three- 
frequency forcing has been used too 
extent. 
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A primary focus of the past 10 years of investigation 
of one- and two-frequency forced Faraday waves has been 
on the role of resonant triad interactions - the lowest or- 
der nonlinear interactions - in pattern selection. A typ- 
ical resonant triad consists of two critical modes and a 
third, damped mode with a distinct wavelength that de- 
termines the angle rcs between the critical wave vectors; 
see Fig. ^ It was originally thought that the damped 
mode should draw energy from the excited modes, creat- 
ing an "anti-selection" mechanism that suppresses such 
triads in favor of patterns that avoid the resonant angle 
[l2j |. However, it has since been shown that the opposite 
effect may also occur. For instance, demonstrated 
for two-frequency forced Faraday waves that a damped 
mode oscillating at the difference of the two forcing fre- 
quencies is responsible for selecting the superlattice-I pat- 
tern reported in 




FIG. 1: Fourier space diagram of resonant triads composed 
of two critical modes (|ki| — | | = fci) and a damped mode 
(k 3 = ki+k 2 , |k 3 | = fe 3 ) with (a) fc 3 < fci, (b) fci < fc 3 < 2ki. 
The resonant angle is given by cos(# res /2) = k 3 /(2ki). 



The richness of the Faraday system is due in large part 
to the vastness of the control parameter space which, 
in principle, is infinite-dimensional. By expanding an 
arbitrary periodic forcing function in a Fourier series 



F(t) = f m e 1 ' 
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it is apparent that one has the freedom to choose the com- 
mensurate frequencies (mw, nu>,pu>, . . .), the amplitudes 
(l/m|,|/rj|,|/p|)---)j and the phases (<f> m , <j> n , tj> p , . . .), 
where <p u = arg(/„). Previous experimental and the- 
oretical investigations have largely been exploratory in 
nature, i.e., they have focused on describing the patterns 
which form at particular locations in parameter space. 
In this Letter, we take a prescriptive approach to pat- 
tern formation and demonstrate how the huge parame- 
ter space may be exploited to control resonant triads. In 
particular, we can either enhance or suppress triad in- 
teractions, largely at our discretion, by an appropriate 
choice of forcing frequencies, amplitudes, and phases. 

Our analysis follows that of [llj which uses the (bro- 
ken) symmetries of time translation, time reversal, and 
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Hamiltonian structure, and is valid near onset in weakly- 
damped systems. Our main result is a table summariz- 
ing, for up to three forcing frequencies, which damped 
modes are likely to be important, the manner in which 
their coupling with critical modes depends on the forc- 
ing parameters and, in many cases, the overall qualita- 
tive effect (enhancing or suppressing) they have on as- 
sociated patterns. These results are then applied to two 
specific examples relevant to experiments. In the first, we 
show how the superlattice pattern selection mechanism 
discussed in 0, 0] may be dramatically enhanced by 
the addition of a third frequency component with appro- 
priate phase. In the second, we offer an explanation for 
why a distorted 8-fold quasipattern observed with two- 
frequency forcing became much more regular and robust 
when a third forcing frequency component was applied 
in the experiments of 0. 

We consider resonant triads composed of two critical 
modes and a third linearly damped mode (which may, 
nevertheless, be forced) oscillating with dominant fre- 
quency SI > 0. The particular values of SI most crucial 
to pattern selection will be determined in the course of 
our calculation. Time is rescaled so that the common 
frequency is one, and we take the m component of the 
forcing to drive the critical modes, which therefore have 
dominant frequency m/2. All other modes are linearly 
damped. We begin by expanding the fluid surface height 
h in terms of six traveling wave (TW) modes: 
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h(x,t) = ^^Z±(i)e l(fe ^ ± ^ t) +c.c. + ..., (2) 
j=i ± 

where w\ — w 2 — m/2 and uj^, = SI. The TW amplitude 
equations respect the spatial symmetries of the problem 
(translation, reflection through k 3 , inversion through the 
origin) and the temporal symmetries 

T T : Zf -> e ±l ^ T Z±, /„ -> e mr f u , (3a) 
n-.Zf^Zf, (t, 7 )--(t, 7 ), /„-/„. (3b) 

Here u denotes any of the frequencies {m, n,p, ■ ■ ■ } and 
7 = 2vk 2 is a dimensionless damping parameter formed 
from the kinematic viscosity v and the characteristic 
wave number fc(m/2), determined by the hydrodynamic 
dispersion relation. Note that T T and k are not symme- 
tries in the usual sense; time translation T T is broken by 
the forcing and time reversal k by the damping. Nonethe- 
less, they are valid parameter symmetries of the problem. 
The equivariant TW equations take the form 

Z+ = vZ+ - Xf rn Z- + Qx{Z±,Z±) + • ■ • , (4a) 
Zt = QZ+ - fiF 2 n Z 3 " + Q 3 {Zf, Zf) + • • • , (4b) 

with the remaining four equations related by symmetry. 
The parametric forcing term F 2 q represents products of 



the /„ whose frequencies sum to 2f2 (i^n = fin when 2S1 
forcing is present). The resonant terms 

fii - (QiZ+ + Q 2 Z^)Z+ + {Q 3 Z+ + QiZz)Z^, 
Q 3 = Q 5 Z+Z+ + Q 6 (Z+Z 2 - + Z^Z+) + Q 7 Z^Z 2 , 

couple the two modes at quadratic order. The res- 
onant coefficients Qe, according to T T , must trans- 
form as {Qx,Q 5 ) ~ e *»~ fi )-, (Q 2 ,Q 7 ) ~ e *(">+n)r 
{QziQiiQs) ~ e~ J ° r . The dependence of the TV- 
invariant coefficients v, A, etc., on 7 and the f u , is also 
determined by symmetry. In particular, k forces most of 
them to be purely imaginary at leading order. However, 
if fci and fc 3 are defined by the local minima of the neu- 
tral stability curves then the imaginary parts of v and g 
(i.e., the detunings) vanish at 7 — 0. We will need only 
the leading term in each coefficient so we write 

v = -v r j, Q=-Qr1, X = iXi, (i = i[u, Qe = iqeFi, 

where v r ,g r > (they correspond to damping terms), 
Xi,fii > (see 14]), qe e M and Fu denotes an appro- 
priate product of the /„ (or unity) . The critical value of 
|/ m | is 0(7) and all \ f u \ are assumed to be related by a 
finite 7-independent ratio, so that ~ 7. 

For a given damped mode, the task is to determine how 
the resonant coefficients Qg depend on the /„, then ap- 
ply a standard reduction procedure at the bifurcation to 
standing waves (SW). The resulting two equations (trun- 
cated at cubic order) describe SW dynamics near onset: 

A x =XA 1 +A 1 (a\A 1 \ 2 + (b + b rcs )\A 2 \ 2 ). (5) 

The equation for A 2 is obtained by switching labels. The 
self-interaction coefficient a and the "nonresonant" part 
60 of the cross-coupling coefficient are C( 7 ). 

Of primary interest here is the contribution 6 rcs makes 
to the cross-coupling coefficient in © as a result of the 
slaved modes Zf . Loosely speaking, if b Tes > the sta- 
bility of patterns involving critical modes separated by 
the angle 9 rcs will be enhanced, whereas 6 ros < has a 
suppressing effect Since we are interested only in 

damped modes that have a significant impact on 0, we 
consider those cases where 6 ros is 0(7) or larger |l7j ]. 

There are then two possibilities for obtaining large 6 res . 
If some of the Qi are 0(1) then 6 res is C(7 _1 ) and domi- 
nates in the limit 7 — > 0; this happens only for = m and 
is in essence a single-frequency phenomenon (the "first 
harmonic" resonance). If SI 7^ m and some of the Qe 
are 0(7) then 6 res is 0(7), and hence comparable to 60. 
This occurs when SI £ {2m,rt,m ± n,rt — m} for some 
frequency n. However, the SI = 2m mode may be ig- 
nored because for Faraday waves it has a wavenumbcr 
k 3 > 2kx and so it cannot satisfy the necessary spatial 
resonance condition. Each of the remaining conditions 
on SI generates a particular type of coupling. For exam- 
ple, with SI = n — m we have F 2 — Fj — f n . If more than 
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one condition is satisfied, as when fl = n — m = m — p, 
there are more coupling terms (here one would also have 
Fi = F 5 — f p ). The maximum number of conditions that 
f2 can satisfy simultaneously is three. 

In addition to the issue of coupling terms, it matters 
whether or not the damped mode is parametrically forced 
at 0(7). This forcing is present when there is a fre- 
quency 20, in F(t) and magnifies the resonance effect as 
it brings the damped mode closer to criticality. In ad- 
dition, it leads to interesting phase dependence as the 
"usual" phase, dictated by /2a, competes with the phase 
of the nonlinear forcing Q3 (see Eqs.0J. 

In Table [I] we give the leading contribution to 6 rcs of 
the important damped modes for different choices of F(t) 
containing up to three frequencies (m,n,p). To simplify 
the expressions therein, we define 

Q'i=9i95, a 2 = 9297, a 3 = 2q 6 (q 3 - q 4 ) , 

&4 = <M7 - 9295, «5 = £\(2qiq6 + 95(93 - 94)), 

a 6 = e A (2q 2 96 - 97(93 - 94)), £a = sign(Aj), (6) 
and the functions 

p (ib \ = (lei + l-fen I sin $) 
201 ' (|<?| a -|rt/an| 2 ) ' 

fi 2 n(^^ 2 )= (|g|si ^rr l t Ql r $2) 

(\Q\ 2 - \Vihn\ ) 

There are four groupings in the table. The first shows 
the five important damped modes and their contribution 
to 6 res when there is only one type of coupling at 0(7) or 
lower and no parametric forcing /2n- In these cases there 
is no (leading order) dependence on the forcing phases. 
In the second section the same damped modes have been 
parametrically forced. The factor X/\q\ is then replaced 
by P2Q,{^) of 0, a strictly positive oscillatory function 
(\q\ > \liifm\ for damped modes) with extrema at $ = 
±90°. Entries in the third section display two types of 
coupling, while the final two cases in the table have /2fi 
forcing as well. Note that equivalent cases can be trivially 
generated from those in Table[l]by switching n and p and 
relabeling, for example (to, n,n — to), = p is equivalent 
to (m, n, to + n), fi = n. 

Hamiltonian structure in the undamped problem has 
important consequences for the results in Table Q] If one 
assumes that at 7 = Eqs. derive from a Hamiltonian 
H through dZf /dt = TidH/dZf (see, e.g., |l|]) then 

91 = 95, 92 = 97, and 93 = 94 = 96- When we allow 
for simple rescalings: (Z^,Z^) — * rj(Zf , Z^ 1 ), {Zf) — > 
£,{Zf) with ij^eK, these conditions relax to qi = rq$, 

92 = ^97, and 93 = 94 = rqg (for some r > 0), implying 



(7) 
(8) 



a\ > 0, a.2 > 0, a 3 = 0, 0:4 = 0. 



(9) 



Relations (0 mean that for simple couplings (the first 
two sections of Table[|J the sign of 6 res is determined, and 
thus one knows if the resonant triad has an enhancing or 



suppressing effect on patterns involving 9 rcs . The O = m 
and Q = m =fc n modes are suppressing, the Q = n mode 
is inconsequential, and the Q = n — m mode is enhancing. 

We now pursue the implications of Table [I] by exam- 
ining two cases relevant to recent experiments. We com- 
pare our theoretical predictions with coefficients calcu- 
lated numerically from the Zhang- Vinals Faraday wave 
equations ^3], which describe weakly-damped fluids in 
deep containers. The calculation (see for details) 
gives us the cross-coupling coefficient b = bo + b ICS at 
the SW bifurcation and is independent of the symmetry 
arguments used here. 

As our first example we consider the superlattice-I pat- 
tern observed with two-frequency (to, n) = (6, 7) forcing 
in the experiments of 0. The wave vectors making up 
the pattern lie on the vertices of two hexagons, one ro- 
tated by an angle 6*h < 30° with respect to the other (see 
Fig. 3 in 01)- ^ was shown in ^(J that the experimen- 
tally observed angle of 6*h — 22° is related to a resonant 
triad at 6 ~ 158° (= 180° - 22°) involving the fl = n - m 
mode. This is the most interesting of the damped modes 
because it gives 6 rcs > and thus acts as a selection 
mechanism. Table [I] provides a means by which to en- 
hance this effect, namely by parametrically forcing the 
damped mode via (to, n, 2n — 2m) = (6, 7, 2) forcing. In 
Fig. Gt. we show b(8) for the phase <& = 90° which opti- 
mizes the stabilizing effect at 9h — 22°. In Fig. ^> we 
show that parametrically forcing the difference frequency 
mode with this phase quadruples the stabilizing impact 
with respect to the two-frequency case used in the ex- 
periments. Note that if the wrong phase ($ = —90°) is 
chosen, the effect of the resonance will actually be dimin- 
ished compared to the two-frequency case. Fig.|2t shows 
the sinusoidal dependence of 6 1CS on $, in excellent agree- 
ment with our predictions. 
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FIG. 2: Effect of $ on the stabilizing = n — m resonance 
that selects superlattice patterns with (m,n,p) = (6, 7, 2) 
forcing, (a) b(9) with the optimal phase $ = 90°; the singular 
region heralding hexagons at 9 = 60° is removed. The large 
dip near 9 = 67° degrees is due to the strongly suppressing 
Q, = m resonance in the first section of Table (b) Close-up 
of b{9) near 9 = 22° with $ = 90° and $ = -90°; the two- 
frequency result (dotted line) with f p — is also shown, (c) 
Spike magnitude versus $ (see Eq. [7J . For these calculations 
|/n|/|/m| =0.75, \f P \/\fm\ =0.1. The fluid parameters in the 
Zhang- Vinals equations of [l2Tl are 7 = 0.08, Go = 1.5. 

As a second example, we consider recent experimental 
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TABLE I: Leading resonant contribution 6 ros in JjsJ for the most important damped modes under appropriate choice of three- 
frequency forcing. Here, m,n,p,S~l > and x € Z + . Each expression for (m,n,p), given f2, is excluded from those of entries 
further down the table. Dots in the first column indicate an arbitrary commensurate frequency. For * the ± follows sign(m — n). 
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results on quasipatterns. It was reported in 7] that 8-fold 
quasipatterns, which were distorted and difficult to ob- 
serve with (m, n) = (3, 2) forcing, became dramatically 
cleaner and more robust with (m,n,p) — (3,2,4) forc- 
ing. An explanation for this observation is provided by 
Table |U Specifically, we find that the = 1 mode forms 
a resonant triad with the critical modes with associated 
angle Tes rj 41°. This is the anglepresent in the dis- 
torted (3,2)-forced quasipatterns in Q- Table [I] indicates 
that with (m,n,p) = (3,2,4) forcing there is a positive 

I/4 1 2 contribution to 6 10S which has a stabilizing effect. 
In our numerical investigations, we find that the stabiliz- 
ing spike in b(6) becomes broader with increasing 7, and 
thus it is reasonable to expect that the stabilization could 
extend to the 45° angle associated with the perfect 8-fold 
quasipattern. Unfortunately, the experimental value of 7 
appears to be too large for a more quantitative compar- 
ison with our theory. 

The theoretical framework developed in this Letter 
should provide some welcome guidance for navigating the 
huge parameter space of multi-frequency forced Faraday 
waves. Table [I] provides a comprehensive list of the im- 
portant damped modes and their effect on pattern for- 
mation (additional possibilities with more than three fre- 
quencies will be considered elsewhere 0). In general, 
the influence of these damped modes, when parametri- 
cally forced, depends on T T -invariant combination(s) of 
forcing phases. Using the "proper" phase greatly en- 
hances resonance effects while the "wrong" phase can ac- 
tually reduce them. We hope that experimentalists will 
use Table|Jas a tool to design forcing functions conducive 
to particular patterns. For example, one might stabilize 
a quasipattern composed of 2N, N > 3, equally spaced 
critical modes by arranging for an f2 = n — m mode with 
wavenumber = 2k\ cos (^^-90°). In the inviscid limit 
an appropriate choice of m and n can be estimated from 
this geometric condition and the dispersion relation 0] . 
Adding a third 2n — 2m frequency would then allow one 



to manipulate the effect by varying <!>. Finally, we em- 
phasize that the scheme used to obtain TableHJis based on 
symmetry considerations and is therefore quite general, 
requiring only that the damping be small. 
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